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A self-consistent formulation is proposed to generalize the HF scheme with the incorporation of 

^»T->~, screening effects. For this purpose in a first step, an energy functional is defined by the mean 

value for the full Hamiltonian, not in a Slater determinant state, but in the result of the adiabatic 

connection of Coulomb plus the nuclear (jellium charge) in the Slater determinant. Afterwards, 

fsj , the energy functional defining the screening approximation is defined in a diagrammatic way, by 

imposing a special "screening" restriction on the contractions retained in the Wick expansion. The 

generalized self-consisting set of equations for the one particle orbitals are written by imposing the 

extremum conditions. The scheme is applied to the homogeneous electron gas. After simplifying 

^—\ ■ the discussion by assuming the screening as static and that the mean distance between electrons is 

«ys ' close to the Bohr radius, the equations for the electron spectrum and the static screening properties 

are solved by iterations. The self-consistent results for the self-energies dispersion does not show 

the vanishing density of states at the Fermi level predicted by the HF self-energy spectrum. In 

this extreme non retarded approximation, both, the direct and the exchange potentials are strongly 

screened, and the energy is higher that the one given by the usual HF scheme. However, the inclusion 

of the retardation in the exact solution and the sum rules associated to the dielectric response of 

the problem, can lead to energy lowering. These effects will be considered in the extension of the 

CZ3 ' work. 

1-i ', 
CS . PACS numbers: 71.10.Fd,71.15.Mb,71.27.+a,71.30.+h,74.20.-z,74.25.Ha,74.25.Jb, 74.72.-h 

0: 

I. INTRODUCTION 

a. 
o ■ 

^j ■ The development of band structure calculation procedures is a theme to which an intense research activity has 

been devoted in modern Solid State Physics. This area of research has a lone history due to existence of important 
unsolved and relevant questions concerning the structure of solids (lH27l. I29J . |30L |43| . One central open problem is the 
connection between the so called first principles schemes with the Mott phenomenological approach. This important 

^^ . approach, furnishing the descriptions for a wide class of band structures of solids, is not naturally explained by the 
assumed fundamental first principles approaches [3|, [a, 0, [jjl, [3(| ■ A class of materials which had been in the central 

VO , place in the existing debate between these two conceptions for band structure calculations are the transition metal 
oxides (TMO) 0, Qjj, GJ, H El HI EJ El • A particular compound which is closely related with the TMO is the 

£T) [ superconductor material La^CuOi in which the first band structure calculations predicted metal and paramagnetic 
characters, which are largely at variance with the known isolator and antiferromagnetic nature of the material [29| . 
Actually, these two qualities are fully believed to be purely strong correlation properties which are not derivable from 
HF independent particle descriptions [lg, Eabij. However, in Refs. [35], |36| a single band Hartree-Fock (HF) study 
was able to produce these assumed strong correlation properties of La^CuO^, as independent particle ones, arising 
from a combination of crystal symmetry breaking in and an entanglement of the spin and spacial structure of the 
single particle states. Therefore, this result opens the interesting possibility of clarifying the connection between the 
Mott and Slater viewpoints in the band theory of the TMO [35], ]36| . It should be underlined that a former HF 
study of La,2CuOi in Ref. [37! ],, also indicated an insulator an AF character of the material. However, the extremely 
large gap (17 eV) which was predicted by that work had a radical discrepancy with the experimentally estimated 
value of 2 eV. However, it could happens that the results of Ref. [37j could show a close link with the conclusions 
of Refs. [33, |36|. Results in the literature supporting this possibility are the ones presented in Ref. [38(, where the 
same recourse employed in Refs. [35], [36| of allowing a crystal symmetry breaking by doubling the unit cell of the 
lattice was employed. In that work, the spin was treated in a phenomenological way by including a spin dependent 
functional in the variational evaluation done. These two procedures allowed the authors to predict an insulator and 
AF structure for the material with a gap of 2 eV, well matching the experimental estimate. 

The results in Ref. [3a, |36[ were obtained by starting form a Tight Binding (T.B.) model with a square lattice 
size equal to the distance between Cu atoms in the CuO planes, by assuming a simple crystalline extended quadratic 
(near the Cu atoms) T.B. potential. The strength of the potential was fixed to reduce the width of the states to be 
smaller than the chosen lattice constant to obey the T.B. condition. Afterwards the metal paramagnetic HF bands 
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obtained when considering the Coulomb interaction in the model, was fixed to approximately reproduce the single 
half filled band crossing the Fermi level, obtained in the evaluation of the band structure of La-iCuO^ in Ref. |29| . 
A curious point was that the value of the dielectric constant that was required to make equal the two bandwidths 
was nearly e = 10 . Therefore, it come to the mind that the introduction of a dielectric constant in the HF equations 
of Ref. [23], might be produced a band gap of nearly 1.7 eV for La^CuO^ being very close to the measured one 
of 2 eV. This indicates that if a HF method could be extended to predict the screening properties, it could has the 
chance of explaining the strongly correlation properties of LaiCuO^ and with luck, of its close related parents, the 
transition metal oxides. The possibility exist that such an approach could furnish a simple way of catching the strong 
correlation properties of these materials giving in this way a clear explanation for the Mott properties, by showing its 
connection with the crystal symmetry breaking and spin-space entangled nature of single particle states [35|, |36| • 

Therefore, in this paper we intend to start generalizing the HF scheme to introduce screening effects. A first 
possibility is considered here by constructing an energy functional E to be optimized by a set of single particle 
states. The functional E for the HF scheme is simply taken as the mean value of the exact Hamiltonian in the Slater 
determinant of formed by the single particle functions over which E is optimized. 

Since we are interested in introducing Coulomb interaction effects, a first general possibility consists in modifying 
the variational Slater determinant of the single particle functions by connecting in it the Coulomb interaction between 
the electrons and with nuclear "jellium" potential by employing the Gell-Mann-Low procedure |40l443| : This definition 
will have variational character, that is, will be equal to a mean value of the total Hamiltonian in one state, and thus 
the result after finding the extremes over the one particle states, should be greater than the ground state energy of 
the system. However, the complications in solving the extremum equations following from introducing the Coulomb 
and nuclear interaction in an exact form will be enormous. However, we are only interested in catching improvement 
on the HF scheme produced by more simple screening processes. Then, it becomes possible to approximate the above 
mentioned state (the one resulting of the connection of the electronic interactions) by restricting the sum over all 
the Feynman graphs to a set reflecting screening properties. Then, the simpler state will be chosen by only retaining 
in the sum given by the Wick expansion, only the graphs for which: a) for any line arriving to a point of a vertex 
and coming from a another different point, always exists another different line that joins the same two points. This 
definition basically retains the diagrams in which are essentially the "tree" diagrams (having now closed loops) in 
which arbitrary number of insertions of the polarization loop are included. The functional E defined as the mean value 
of the Hamiltonian is also an upper bound for the ground state energy. However, this functional expression although 
simpler, yet will lead to slightly cumbersome procedure. Thus, we defer its study to eventual future discussion. In 
order to further simplify the form of the functional, but yet retaining screening effects, we even more simplify the 
functional by imposing the same screening approximation to the whole sum of Feynman diagrams which defines the 
functional E in which the Coulomb and nuclear interactions are adiabatically connected. This definition leads to 
scheme which is formulated in very close terms with the usual discussions of the screening approximation in Many 
Body theory [H. 

Up to now, we had not been able to show that this approximation also has a variational character. However, at 
least it follows that the Feynman diagrams defining the scheme represent allowed physical process, which suggests 
that the procedure can improve the HF results. If the variational character turns to be valid, the result for energy 
will also furnish an upper bound for the ground state energy. This issue will be considered elsewhere. 

After defining E, the generalized HF equations are derived for the general case of inhomogeneous electron systems. 
They include exchange and direct like terms, but also show additional contributions being related with effect of 
the Coulomb interaction and nuclear potential on the kinetic energy of the electron which vanish for homogeneous 
systems. The equations are also written for the case homogeneous electron gas in which the full translation symmetry 
simplifies them. In this case, as in the usual HF scheme, only the kinetic and exchange terms are non vanishing. In 
this first exploration, it is assumed that the screening is static and the density is chosen to be one electron per sphere 
of radial dimension equal to the Bohr radius. Under these conditions, the self-consistent equations for the electron 
spectrum and the static screening properties are solved by iterations. The result for the dispersion does not present 
the vanishing density of states at the Fermi energy appearing in the fist in the HF self-energy spectrum [44H48J . 

The work proceeds as follows. In Section II, the generalization of the HF functional and equations of motions for 
the single particle orbitals are presented. Section III then considers the application to the homogenous electron gas. 
In the conclusions the results are reviewed and possible extensions of the work are commented. 



II. GENERALIZATION OF THE HF SCHEME TO INCLUDE SCREENING EFFECTS 

A. Hartree Fock approximation 

As it is known the HF approximation is obtained after minimizing the mean value of the full Hamiltonian of the 
system H in a Slater determinant state | <I>o) formed with the elements of a basis of wavefunctions (fk- For future 
purposes, let us define a free Hamiltonian Hq in terms for these functions in the form 

H =^e k a+a k , (1) 

k 

where a~£ and a k are the creation and annihilation operators of the functions y>fc(x) and the energy eigenstates {e k } 
are to be defined by the optimization process of the energy functional 

E[<p] = ($ |# ^o> -$>«(/ dx<4(x) W (x)-<^), ( 2 ) 

kl J 

(x 1 ,x 2 ,...,Xjv p |$o) = —^ Det ifk(^i)} (3) 

where N p is the number of electrons. The bold letters will represent the combination of the position x and spin 
s = ±1 coordinates asx= (x, s) and the integrals will mean 
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s=±l' 

In addition, in order to simplify the expressions when time dependence will be considered we will define also space-time 
coordinates x = (x, t) = (x, s, t). 

The system to be considered will be a set of electrons which move in the presence of a nuclear potential V nuc 
and interact among themselves with the Coulomb two body potential v. After adding the kinetic energy T of the 
electrons, the total Hamiltonian is given in the Schrodinger coordinate representation as 

H = T+V nuc + v, (4) 

Np fi 2 8 

Np f -+ e 2 

Vnuc = -^ / d ~f Pnuc{~£) ,-A. _ -^ , , (6) 

^ 1 e 2 

and, in what follows, in order to simplify the notation, all the operators acting in coordinate space appearing, will be 
supposed to have a spin structure, which when not written, are assumed to be the identity matrix in the representation 
of the functions on which the operators act. Imposing the extremum conditions on E over the variations of the 
functions, the HF equations for determining the single particle states follow. The fact that, on the extremum values 
of the single particle functions, it is satisfied the E = ($ | H ^o) 1 implies 

E > E g , (8) 

where E g is the ground state energy, which for a Hamiltonian system gives the lowest mean value possible in any 
state. 

Let us now intend to construct a generalization of the HF method including screening effects. For this purpose, 
initially, let us define the functional to be optimized in an alternative form assuring that the result includes Coulomb 
interaction effects. A possible form of this kind is the following one 

E = {M{U c a (0, -oo))+ H E£(0, -co) ^o) 



Po|(tfS(0,-oo))+tfS(0,-oo)|&o) 

fe |t^(oo,0)ff E/°(0,-oo) po) 

p \UZ{oo,0)Uz(0,-oo)po) ' 
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in which [/£ is the evolution operator connecting the Coulomb interaction among the electrons plus the nuclear 
potential V nuc + w, in an initial previously defined Slater determinant | $q) (the ground state of the free Hamiltonian 
Ho) at a large time in the past 

U c a (0,-oo) = J2^T dtl dh - dt n T[H c a ( tl ),H c a (t 2 ),H c a (t 3 )...,H c a (t n )i 

_n J —00 J — 00 J — 00 

= T[eM~J_ dtH c a (t))}, (fO) 

H c a (t) = exp(lff i) H c cxp(-^H t) e^~ a I * l\ (ff) 

# c = K„ c + v, (12) 

where iJ c is the sum of the Coulomb and nuclear electron interactions in a time independent Schrodinger repre- 
sentation, a is the Gell-Mann-Low small parameter which implements the adiabatic connection of the interaction 
when taken in the limit a — > + . The a dependent exponential factor a p pea ring in the definition of the interaction 
representation operators H^(t) above, tends to one when a vanishes [40l. |4U |43| . The T symbol appearing is the 
time ordering operation which orders the operators from right to left for increasing time arguments |40l . l4l[ |43| . 
For the Hermitian conjugate operator U^ + it follows 



(E£(0,-oo)) + = J2^~ dtx dt 2 ... / dt n T* [H c a {t x ),H c a {t 2 ),Hl{H)...,H c a {t n )] 



n=0 

= .T*[eM^fjtH c a (t))} 

if if 00 

= T[exp(-y dt*H c a (t*)))=T[eM-j:j dt*H c a {t*))] 

= J7=(oo, 0). (13) 

Noticing that H is equal to its interaction representation Hit) = exp(jrHot) Hexp(—iHot) e^ a I * D at t = 

ff (0) = T(0) + 14 uc (0) + v{0), (14) 

it follows that in formula ((9]), the operators are temporally ordered. Therefore, the Wick theorem can be applied. 
The sum of the expressions associated to all the Feynman diagrams giving the denominator, coincides with the sum 
of all the possible terms associated to disconnected (unlinked) graphs of the numerator, multiplying the contribution 
of a given arbitrary connected graph. Thus, this common factor cancels [40, [4j], |43| . The expression for E in this 
general form is given by the sum of the contributions being associated to all the diagrams being connected through 
their vertices and lines to the vertices defined by the operator H. In this case, since E is again a mean value of H in 
a given state, it follows that E should be greater than the ground state energy E g 

In spite of this property, the complexity of the sum over all the connected part will surely make the optimization 
problem an impossible to be solved one. Thus, it is needed to further simplifying the functional incorporating 

the screening effects. 

B. Screening approximation 

As it was just mentioned, the general form of the variational problem defined by @ is too complicated for to be 
usable in concrete evaluations. Let us consider in this section two simplifications. They will be called as "screening" 
ones, and will be implemented through an approximation of the Wick theorem. In the Wick expansion, the time 
ordered products of creation and annihilation operators are expressed as the sum of all normally ordered products 
including an arbitrary number of contractions between any pair of operators entering in the product. The terms in 
this ex pan sion are related by a one-to-one mapping with all the diagrams generated by the Feynman rules of the 
theory g3,0,lil • 

The first "screening" approximation to be considered will be one in which a variational state | $ s ) = U& \ $0} 
is defined after applying the Wick expansion to the operator E/£(0, — 00) entering in the definition of the state 
t/„(0, — 00) p>o)- In the Feynman diagram expansion of the operator only the graphs will retained in which, for 



any particular line joining two different ending points of any vertices, another different line always exists connecting 
the same two points 4l| . We will describe this reduction of the diagrams associated to a time ordered operators A 
by writing [AJiy,s- Then, the variational state | $ s } will have the form 
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and its Hermitian conjugate state satisfies 
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It is clear that since the E functional is defined as a mean value in the given normalized state, the minimization of 
E process will be variational, that is E > Eg. Since the criterion for eliminating a contribution is only determined by 
the contractions between two linked vertices, the application to the screening approximation to a product of two to 
disconnected graphs G\ and G 2 in which all the operators are contracted, gives the product of the two independently 
approximated graphs. That is, the same property as for the Wick expansion 



[G\G2]ws — [G\]ws[G2] 
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The approximation done in defining the state eliminates all the graphs in the Wick expansion of the state 
t/„(0, — 00) |$o)- But when the scalar product defining E is taken, a new application of the Wick expansion to 
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the full operator Ua HU a leads to diagrams not satisfying the defined screening approximation. That is, diagrams 
exist for which not all the internal fermion lines appears in polarization loops and tadpoles. This fact is not an 
obstacle of principle, but it leads to an structure of the connected to H sum of Feynman diagrams in which somewhat 
complicated three and four loops diagrams appear. Then, in this work will not use the approximation in this form. 
Rather, in order to retaining contributions consistently being given by diagrams showing the imposed "screening" 
graphical restriction, the functional E will be defined in the form 



r(c) + 
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E = 
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in which, the modified Wick expansion [ ] w $ is applied to the whole operators Ua HU a and U a 
After expanding the exponential defining the operators, the numerator has the expression 
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Hx 

TK(*i),^(4),^(*3)-^a(C)]]vr5^0), 

and the denominator has a similar one, obtained by omitting the H operator. 

The figure ((T|) is a pictorial representation of the terms associated to the product of the power (m 
expansion of U a and the power (m + m^) of the expansion of the operator U a in formula p^|) for the E functional. 
Each square represents one interaction operator v + V nuc evaluated a given time, which grows from right to left in the 
picture. The illustrated term represents a particular contribution to the approximated Wick expansion adopted. In it, 
the white squares represent vertices which have contractions between the operators defining them. Such contractions 
are in such a form that the graph associated to these vertices and contractions are not connected to the vertices in the 



(19) 



-mi) in the series 



D...DD.. .dHdD ...DD-D 

n 2 1 m 2 2 1 12 m, 1 n, 

t>Q t=0 t<0 



FIG. 1: The figure illustrates the terms associated to the product of the power (m + mi) in the series expansion of Ua and 
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the power (ni + mz) of the expansion of the operator U a in formula (|19|l for the E functional 



operator H. The white channels joining such vertices represent the set of contractions in the illustrated contribution. 
Alternatively, the black squares represent all the vertices of the interactions v + V nuc being connected among them 
and with vertices in H by contractions. The square labeled with the letter H represents any of the vertices in the 

(c) 

total Hamiltonian. The number mi of black squares come from the expansion of U a and the number rri2 from the 

(c) + 

expansion of U a . The black channels in this case symbolize the set of contractions between operators associated to 
these vertices. The fact that all the terms in the Wick expansion are normal orderings including contractions, allows 
to reorder in the time direction all the square boxes in an arbitrary way. This follows because each vertex has an even 
number of these operators and the normal ordering are even under a permutations of groups of even numbers of such 
operators. Therefore, the same contribution to the result of the expansion will be obtained from the 



(wi+mi)! 



terms 



appearing the expansion of the exponential operator U a , each of which has a coefficient , . .,, . In the same way, the 



(ni+mi). num ] : , er f identical contributions in the expansion of U a changes the original coefficient 
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Thus, each connected contribution will be multiplied by the sum of the same contributions defining the denominator 
($o|[ka U a jwsfaa)- in which H is not appearing. Therefore, the functional E is defined by the sum of all the 
contributions associated to the Feynman diagrams only including "screening" contributions being connected to the 
operator H. It is not ruled out this second type of screening approximation also gives a value of E being an upper 
bound for the ground state energy. However, we had yet found a proof of the property, and it will be expected to be 
considered elsewhere. In spite of this fact, since the HF approximation furnishes good results for the energy of many 
systems of interest, leads to the expectation that the inclusion of the screening effect can provide improvements of the 
HF results for the total energy as well as for the excitation spectrum at least in some situations of physical interest. 



C. Analytic expression of the functional 



Let us write in more detail the Hamiltonian of the system as written in second quantization, and the rules for the 
Feynman diagrams in the considered general many electron problem. The definition and notation follow the ones in 



Ref. j42[. The Hamiltonian in second quantization has the form 

H(t)=T(t) + V nuc (t) + v(t), (20) 

T(t) = -Jdx^ + (x)d?*(x), d i = -^-, (21) 
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V nuc {t) = - / d -f Pnuc C?) -j e __g I d x *+(s)*(a;), (22) 

«(*)= /da /riz'V(^,8,t)\l/ + (t,a,t) ^"J $ + (t,g,t') ^+C^,a,t), (23) 

*(x) = y^Qfc(t) ^fe(x) = ^]a fc exp(— -efct)(/?fc(x), (24) 

*+(*) = XX (*) ¥>fe( x ) = E a fe exp(le fe t)^(x), (25) 

fe fe 

The above Hamiltonian leads to a diagram expansion as follows. The directed continuous lines represent the free 
fcrmion propagator (contractions) associated to the Hamiltonian Hq 

, [dw^i <p k (x)<p%(x.')exp(-i w(t-t'-S)) 
U(x,x ) - / — > J — 

J 2tt *-£ (w - Ti e k + i a sgn(e k - e f )) 

= ]T <p k (x)<pl(x!)(0(t -t'- 6)6(e k - e f ) - 9(t' - t + 8)0{e f - e k ), 
k 
5^0+, 

where, as defined before, a is the small parameter introduced for connecting the interaction and 6 is small time 
splitting which is required to well define the time ordering for equal time operators in each interaction Hamiltonian 
according to the convention in Ref. [42j. In particular it assures the appearance of the usual HF terms in the 
discussion in the limit of vanishing polarization effects. The wavy lines with two ending points represent the vertices 
of instantaneous Coulomb interaction potential 

in which a factor (— 4) multiplying the Coulomb potential and a delta function reflect its instantaneous character. A 
circle with a central letter n as a label and with a point attached to a Coulomb wavy line ending at point x (at 
which two fermion lines attach, one coming and another outgoing) represents the one body nuclear potential in (|22p 

Hb/^^^pn^r (27) 

The kinetic energy vertex (|21|) 

" 5> < 28 > 

is represented by a circle with a central letter T joint to a point on which one fermion line arrives (in the sense of its 
attached arrow) and another line departs. The Laplacian operator acts on the coordinate argument of the fermion 
propagator at which the fermion directed line arrives. Finally a black circle will represents the sum of the (negative) 
electron charge density and (positive) charge density of nuclear particles. 

The generic types of connected contributions satisfying the screening approximation are illustrated in figure 2. 
The first series corresponds to the addition of an infinite number of self-energy insertions of the fermion polarization 
loops in the Coulomb interaction line of the standard HF exchange term. The grey circles with attached wavy lines 
represents screened potential resulting from the described insertions. The second series corresponds to the sum of an 
arbitrary number of such insertions of the polarization loop in the Coulomb interaction line of the standard direct 
term in the HF scheme. As defined before, this total density is the electron particle density minus the density of the 
nuclear charged particles. This direct term includes the geometric series of polarization insertions in the Coulomb 
interaction between: the electron charges among themselves, the nuclear "jellium" charges among them and the 
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FIG. 2: The graphic representation of the sum of the Feynman diagrams associated to contributions linked with the vertices 
of the total Hamiltonian. The insertion at the bottom symbolize the total particle density (black disk) as defined by the sum of 
the electron density (the closed propagator line) minus the density of charged particles in the nuclei (the circle with an interior 
letter "n"). For homogeneous systems in which the jellium is homogeneous and the electron full translation symmetry is not 
broken, this total density vanishes. 

interaction between these two kind of charges. The last four terms combined correspond to screening contributions 
to the mean value of the kinetic energy, which are mainly determined by the total charge density. Each contribution 
having a given number I of fermion polarization loops in its diagram has a sign (—1)'. The geometric series over all 
the insertions in the graphs of figure 2 reproduces the usual definition of the screened Coulomb potential [33, |42[ . 
The polarization loops in the graphic have the expressions 



P(x,x') = —G(x,x')G(x',x) 

= J <MxM( x 'Vfc(x'M(x)expH w(t-t'))x 
2 /fc — fi 



h(w + i(e fe - ei) + i a(f k - fi) 
fk = 0(e f -ejfe), 



(29) 
(30) 



and the screened Coulomb potential arising form the summation of the geometric scries of polarization insertions is 
given as 



v s (x, x') — v c (x,x') — / dxidx2V c (x 7 xi)P(xi 7 X2)v c {x2,x')+ 

dxidx2dX3d,X4V c (x, Xl)P{xx,X<i)v c {X2, X3)P(x3,X4)v c (x4, x) + 



(-1)" / dx 1 dx 2 dx 3 dx 4: ...dx2n-ldx2nVc(x, Xl)P(Xl,X2)v c (X2,Xa)P(X3,X4.)v c (X4, x').. 



V c (x,X2n-l)P(X2n-l,X2n)v c {X2n,x') 
1 



(Vc 



1 + Pv, 

1 



■)(x,x') 



v c )(x,x') 



1 + VcP 

= I dx\ e(x,xi)v c (xi,x r ), 



(31) 



in which the usual definition e(x, x ') of the dielectric function in the screened approximation appears [39J. 

The screening effect on the kinetic energy can be expressed in terms of a correction to the electron propagator 
including all the self-energy insertions of the total potential generated by the total charge density pt, as follows. This 



screening modified Green function 

G s (x,x') = G(x,x') + f dx 1 G(x,x 1 ){^-)V t (x 1 )G(x 1 ,x r ) + 

— % — 2 

dx2G(x,xi)(—)Vt(xi)G(xi,X2)(—)Vt(x 2 )G(x2,x') + ... 

dx 2 G(x,xi)( — )Vt(xi)G(xi,x 2 )-( — )Vt(x n )G(x n ,x')... 

= (G -i- )(x,x') 

= ( 1 -^^G)(x,x'), (32) 

V t (x) = / dxv s (x,x)pt(x), (33) 

p t (x) =G(x,x+)+ / o„(x), (34) 

x+ = ("^,t + <5), (5 = 0+. 

where the addition of a small infinitesimal positive number 5 to the time component of a four-vector x has been 
indicated as x + — ( x , t + S), 6 — + . Finally the analytic expression for the functional E can be written in the form 

E = — — / dxdx'G(x 1 x')G{x l 1 x)v s (x 1 x') + 

- I dxdx' p t (x)v s (x, x')p t (x'))+ 

dxdx'(^\7 2 x G(x,x')\ x ,_ >x+ )G s (x' 7 x) 7 (35) 

in which all the participating elements are functionals of the basis functions ifk through the propagator G(x,x') as 
follows 

v s (x,x') = (v c — l —){x lX ') (36) 

1 + Pv c 

G s (x,x') = (G l — )(x, x') (37) 

P(x, x) = -^G(x, x')G(x', x) (38) 

h 

Vt (x) = T i dx v " (x ' x>)pt ^ (39) 

p t (x)=G(x,x + )+p n (x) (40) 

[ dw^ i y fc (x)^*(x / )exp(-z w(t-t' -J)) 

<jr(£,£j — / — > j — (41 J 

J k (w ~ jre k + i a sgn(e k - e/j) 

sgn(a) = 6[a] - d[-a\. (42) 

Let us write the variational equations for the single particle wavefunctions ifk(x) in next subsection. 

D. The equations for the single particle states 

As usual, the set of equations for determining the basis functions (fk{x) will be written as the vanishing derivatives 
of a Lagrange functional with respect to all the functions, by also imposing on them the orthonormality constraints. 
Since E is a functional of the ifk{x) only through its dependence on the fermion Green functions, the chain rule for 
the functional derivatives is helpful. Then, the following expression for the derivative of G with respect to any of the 
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conjugate one particle wavefunctions can be employed 



<Vfc(* 



-Mx,x') = Vk (x) «$( 3 >(x" - x')/ fe (t - t') 



f k (t - if) =(0(t -if - 8)9{e k -e f )- 9(t' - t + 6)9(e f - e fc ): 



.£fc , 



.cxp(-tj(t-t' -5). 



(43) 
(44) 



That is, the functional derivative of G with respect to tpk(x), is proportional to the wavefunction at the same 
quantum number, as multiplied by an energy dependent factor of unit absolute value and by a Dirac Delta function. 
The Delta function when integrated over an internal variable x' will define x" as the external point of the diagram. 
The point x will be the one at the other extreme of fermion line where the arrow arrives to the point to in which 
the wavefunction ip k is evaluated. Employing the Lagrange multipliers method to find the extremal of E by also 
imposing the orthonormality constraints on the functions Lp k , the effective functional Ei, ag to be used for writing 
the equations for the extremals, will be 



ELag =E~y^ y Xkl ( / dx (fl(x)(pi(x) - S k l) 
Ijtk •* 



(45) 



where X k i are the Lagrange multipliers of the imposed orthonormality constraints. Note that, the not yet specified 
energies e k defining the free Hamiltonian _ff , are not yet related with the Lagrange multipliers A^. Then, the 
above rule for the derivative of the propagators, allows to write the extremum condition for EL ag in the form 



8<Pk(* 



SXhW 



E^ag — 



/ dx! H Fock (x,x r )ip k (x.') -^A fe ; ipi(x), 



E L ag = dX (fl(x)(pi(x) - Ski = 0, 



in which the modified Fock kernel including screening effects is written in the form 

H Fock (x, x') = T(x, x') + V c k x (x, x') + V* r (x, x')+ 

T/ fc (x,x') +« V^ v (x,x>) +( 2 ) V*_ v (x,x'). 



(46) 

(47) 



(48) 
(49) 




FIG. 3: The picture illustrates the set of generalized HF equations for the single particle states. They follows after taking 
the derivative of the energy functional over the conjugate of the single particle wavefunctions at the space point represented 
by the external left end of each graphic. The squares with an interior ip represent the single particle orbitals being searched. 
The figures connected with them represent the various kernels defining the generalized Fock operator in expression 1481 The 
dashed lines represent the time and quantum number dependent functions fk(t — t'). The time integrals over these functions, 
in combination with the time invariance symmetry implements the time independence of the Fock operator. 
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The various kernels appearing in the above expression have the explicit forms 

h 2 



T(x,x' 

Vdir ( x , x 

T/ fc (x,x' 
v t (x' 



2m 



e(e/-e fe )V^(x-x'), 



= - / dt di7/c(^-i') G ( x : i ) x '^ / )^(x,t,X4,t4)£(x4,i4,x',i / ) J 
= ~6»(e/ - e fc )v t (x)6(x - x') 

dt dt'f k (t - t')v t (x)G(x, x')v t (x')pt(x'), 



v s (x,x")p t (x")dx 



(50) 

(51) 
(52) 
(53) 

(54) 



V T fc (x,x') = -<5(x-x') / dx'v s {x,x 2 )G{x 2 ,x 3 ){- — )V\G{x 3 ,x 2 )) 

«V^U(x,x') = - J dt J dt'dx 2 dx 3 f k (t-t')v s (x,X2)G(x2,x 3 )(-^)V^ 3 G(x 3 ,x 2 )x 
G(x,x')vt(x')+ 



(2) ^(x,x') 



dt I dt' f k (t -t')(-— )V^G(x,x')v t (x') 



dt I dt' Vt ( x )f k (t-t')G(x,x')(-—)V^ 



(55) 



(56) 



dt I dt'dx 2 dx 3 v t {x)fk{t -t')G(x, x')v s (x',x 2 )G(x 2 ,x 3 )(--— )V^> G(x 3 ,x 2 ) (57) 

2m 3 



(58) 



The equations for the single particle wavefunctions are graphically illustrated in figure [3] It should be noted that 
the electron, nuclear and total particle densities are constant, thus the total potential Vt is also time independent. 

At this point we will specify the spectrum of the free Hamiltonian Ho. The most natural way we estimate for 
choosing it, seems to require that these energies coincide with the solutions for the diagonal components of the matrix 
of eigenvalues Am, which are defining by the extremum problem (as functions of the spectrum {e^}). This selection is 
suggested, by the fact that in this case, the free quantum mechanical problem defined by Hq, will have an energy, at 
the extremum solution for the single particle wavefunctions, which exactly coincides with the value of the optimized 
functional E. 



E. The homogeneous electron systems 



Let us assume now that the electronic system under consideration is homogeneous. Then, the quantum numbers 
k can be taken in form k = p = ( p , <r z ), in terms of the spacial momenta and the projections of the spin in a given 
direction. The normalized wavefunctions in a large box of volume V are the plane waves. The amplitudes i^pwill have 
a spin dependence implicitly assumed 



<y5 P (x) = —tp p exp(i ~f.~£) 



(59) 



Since the Lagrangian multipliers with different indices vanish as follows form the Lagrange extremum conditions, 
after defining A p = A pp = e(p) = k, the HF equations take the simple form 

H F ock (pV p = A p tp p , (60) 

H Foc k(p)=T(p) + V? x (p), 

where the Hamiltonian is the sum of the following constants defined by the Fourier transforms of the potential 

T(p)^%- £p )^, (61) 

Kx(p) = -f I & I dr expH frf) I dt dt7 fc (i-t')G(x,t,x , ,i')« s (x,i,x 4 ,i4)£(x4,i4,x',i'), (62) 



e(x,x') = ( 



1 
l + v c P 



)(x,x'). 



(63) 
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Note that various potential terms in the general inhomogeneous case are not appearing because in the homogeneous 
electron systems, the "jellium" nuclear charge should rigorously cancel the electron charge for the system to be 
homogeneous. The vanishing of one of these terms, in which the total density is not appearing, also follows after 
being directly evaluated. The electron Green function in momentum space is 

nt \ cxp(i w 5) 

G(p,w) = - j -p '- -, (64) 

l(w - j-e p + i a sgn(e p - e/)) 

and the screening potential as a function of the polarization is given by 

v s (p) = v c (p) - v c (p)P(p)v c (p)+ 

l>c(p)P(p)Vc(p)P(p)Vc(p) + ••• 

(-l) n v c (p)P(p)v c (p)P(p)v c (p)... 

Vc(p)P(p)Vc(p) 
_ Vc(p) 

1 + P(p)v c (p) 
= e(p)v c (p), (65) 

47re 2 

Finally, the Fourier transform of the polarization function can be evaluated as follows 

P(~ct, w ) = \ I d (^ - ~£')d(t - t') exp(-i if. (it -lt') + i w.(t - t'))G(x - x')G(x' - x) 

■ k -\ 12 J k J k+~cf 

= ^Vh( W + fa-e t+t )+i a (f t -f t+ - t ) (6?) 



d k 2 J~£ j 



k •' k +-$ 



(2nf h(w + \(^- e t+ ->) + i a (f t f t+ ->) ' 



(68) 



dt 



(2tt) 3 ^ V 

k 



Ef- ( 69 ) 



F. Iterative solution for the electron self-energies in the static limit 

In this subsection we now iteratively solve the set of self-consistent equations for the electron orbitals. In order 
to simplify the evaluation, we will assume the validity of the static approximation. That is, the polarization at 
all frequencies P(q,w) will be substituted by its value at zero frequency P(q,0) In other words, the effective 
potential will act in the same instantaneous way as the Coulomb potential does. In what follows, in order to simplify 
the notation we will not use the combined space-spin momenta notation p—(p ,s) which is however helpful in less 
symmetric problems. 

The self-consistent connection between the quantities can be described as follows. Firstly, the screened Coulomb 
potential can be evaluated from the frequency independent formula 



TTpWfiMT)' 



v s (p) = v.0) - : , nf Tn 7^T > (70) 



in terms of the known Coulomb potential v c and the frequency independent polarization expression 

prf o) = - ( — ° {kf ~ k) ~ ° {kf -^ + ^ (71) 

h J ( 27 3 \{e k e t+t ) + i a (6(k f k) 0(k f |f + f |) ' 

which is fully dehned in terms of the energy spectrum {e^}. Then, this set of energies can be determined again in 
terms of the screened potential v s by means of solving the Fock Hamiltonian equations 
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H F ock(~f)<f^ = e 1 ? <£jh (72) 

HFo ck (t) = ^ + v£(t), (73) 

Vl{t) = J ^0(kf - k)0{k } ~\k+ t\)e(t)v 8 (f). (74) 

In a more detailed form, firstly we choose the Coulomb potential v c as the first approximation for v s in the 
Fock Hamiltonian equations, to determine the first approximation for the energy spectrum. Next, the single particle 
energies defined in this approximation, are used to evaluate the polarization P(p , 0), and with this quantity, a new 
approximation for the screened potential v s is calculated. Afterwards, knowing w s , the same cycle is repeated again 
iteratively up to arriving to stable values of the three quantities P( p , 0), v s {p) and e-^-. 

For the homogeneous electron gas the Fermi momentum kf , defined by 

will fully determine the properties of the solution. It can be defined as usual, the mean radius per particle r$, in 
terms of the volume of the gas V and the number of particles N p by V — Np^-r^ . Then, the Fermi wavevector kf 
expresses in terms of ro in the form 

1 



Given the Bohr radius expression as = — i, the value of kf that was be chosen for the calculation was 

k f a B = 7TTT— = L 

That is, the electron gas density is assumed to be close to one electron within the volume of a sphere having as its 
radius as- The expression of the self-energies in the first step of the iterative process, is obtained assuming that the 
v s coincides with the Coulomb potential, that is v s ( k ) = v c (k). Then, the exchange ^V^(p) and the total eu' 
energies coincide with usual HF results 

{1) V10) = J -^6{k f k)9(k f -\t+ t\)v c (t) (75) 

e 2 k 2 f — p 2 kt + v 

= —^—^ log(l^l) + 2k f ), (76) 

Zir p kf — p 

M.5g!-|l ( ^ l06( |fc±£ (77) 

p 2m 2ir p kf — p 



c 



Afterwards, the first iterative value of the polarization P^ (p, 0), was evaluated with the use of (fTTj) by substituting 
in this expression the spectrum t-J in (1771) . The values of P^\p, 0) obtained were then employed to evaluate the 

first iterative values of the screened potential vi (p) through formula (|TD|) . 

Henceforth, the function vi (~p), allowed to start the cycle again, to find the similar quantities in the second step 
eL. , P( 2 '(p ,0) and v s (p) through the same procedure followed before. 

The described cycles were repeated three times, after which the evaluated quantities closely approached to self- 
consistent values of the functions V e £ ( p ), e-^-, P{ p , 0) and v s (p). The third iteration for the screened potential 
allowed to evaluate one further iteration for the energy, which explains the number of four plotted curves in figure 
SJ The results for these quantities in the various steps are illustrated in figures HJUJ The plots in each of these 
pictures show the values of the quantities in the few steps which were required to attain convergence. It followed 
that the energies, and screening potentials rapidly attain their self consistent values after two steps of the iteration 
process. In all of the pictures 21 [5] and [7J the differences between the second, and third iterations can not be noticed. 

Figure 0] illustrates the outcome for single particle energy spectra for the four iterations done in the solving the 
generalized HF equations. The system considered corresponds to an inter-electron mean distance being close to the 
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1.0 



FIG. 4: The figure shows the single particle energy spectra for the four iterations done in the solving the generalized HF 
equations. The system considered corresponds to an inter-electron mean distance close to the Bohr radius. The energy 
dependence in the first step of the iterative values of the energy coincides with the HF result and corresponds to the lower 
curve. The convergence of the iterative process is rapid, as indicated by the almost coincidence in the spectra for the next tree 
iterations, which are given by the three upper almost identical plots. 



a 0.00105 




0.00090 ^^ 

0.042 0.043 0.044 0.045 0.046 0.047 

P_ 

k, 



FIG. 5: The figure illustrates the small difference between the resulting energy spectra of the second and third iterations done 
in evaluating the energies. The spectrum following in the fourth iteration is also shown, but is not evident, because almost 
coincides with the result of the third iteration. 

Bohr radius. The energy dependence of the first iteration, coincided with the HF result. Figure [5] gives an augmented 
picture of the energies of the last two iterations in the low momentum region. It is shown to evidence the rapid 
convergence of the procedure. The energy associated to the second, third and fourth are shown, but the third and 
the fourth, do not differentiate in the scale of the figure. 

Figures |B] and [7] show the results for the polarization P(p, 0) and the screened potential v s (p). The plots indicate 
that the amount of screening at the considered electron density is high. This can be also noticed from the spectra 
shown in figure 2J which evidences a large reduction of the exchange energy of the usual HF solution. However, since 
the screening properties in the present exploratory study, had been assumed in the static approximation, is not clear 
whether or not this effect will also appears when the exact dynamic screening effects will be considered. This question 
is expected to be considered elsewhere. 

The calculations done in the work referred to homogeneous systems and a very particular value of the density, in 
order to simply illustrate the implementation of the procedure. With this same purpose the strong static approximation 
for the polarization properties was also assumed. However, a further investigation of the homogeneous systems is yet 
a subject of appreciable physical interest and it is expected to be considered elsewhere. In particular, the inclusion of 
the dynamic screening effects could modify the results for the energy, which are higher than the HF ones. Finally, it 
can be noted that the derivative of the energy dispersion at the Fermi momentum kf is finite in the solution, thus 
eliminating this known difficulty of the HF self-energy spectrum |4J| . 



III. CONCLUSIONS 



A generalization of the self-consistent Hartree Fock procedure including screening effects was introduced. Expres- 
sions for the HF like functional and its system of self-consistent equations are presented. For inhomogeneous electron 
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FIG. 6: The figure shows the results for the static polarization P in the three iterations in which it was evaluated. Again, the 
difference between the second an the third iteration (corresponding to the two superposed upper curves) can not be noticed in 
the scale of the graphic. 
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FIG. 7: The figure shows the three iterative results evaluated for the screened potential. As for the other plots the differences 
between the second and third iterations (associated to the coinciding lower two plots) can not be noticed. They indicate that 
the amount of screening predicted for the considered electron density is high. This can be also noticed form the spectra shown 
in figure [4] which shows a large reduction of the exchange energy of the usual HF solution. Since the screening properties had 
been assumed in the static approximation, is not clear whether this reduction of the exchange energy will also appears when 
the exact dynamic screening effects will be considered. 



systems, the equations include additional terms to the standard direct and exchange potentials. The presence of the 
nuclear charge is included in the scheme, seeking to allow its further application in the developing of band calculation 
procedures. To start exploring the implementation of the analysis, it is applied to the homogeneous electron gas. The 
system HF equations is then iteratively solved in the further simplified case of assuming that the dielectric properties 
are taken in the unretarded approximation. The electron density is hxed to be approximately correspond to one Bohr 
radius mean inter-electron distance. The system of equations is solved iteratively by approximating, the screened 
potential by the Coulomb one in evaluating the hrst iteration for the spectrum for the electron energies. Afterwards, 
the screened potential and the new spectrum determined by it, were recursively calculated. The results indicate 
that in the non retarded approximation, both the direct and the exchange potential are strongly screened. Therefore 
the energy spectrum behaves very closely to the free electron one. Thus, the the non retarded approximation does 
not produce an improvement in lowering the energy of the standard HF procedure. The possibility for such an im- 
provement was a central motivation for the present work. However, the non retardation assumption is a dynamical 
simplification which is not necessarily compatible with the (not shown here, but expected to be valid) variational 
character of the scheme. Therefore, the hope exist, that after performing a similar iterative solution, but without 
assuming the static approximation for the fermion loops, can produce an energy competing with the the standard HF 
one. This iterative solution, is not very much difficult than the one done in this work. Thus, if it becomes able to 
furnish an energy competing with the HF one, the method could become a helpful one. The existence of sum rules 
for the frequency dependent dielectric quantities rises the expectation about they could incorporate lowering energy 
contribution, possibly becoming able to contest with the exchange energy in the usual HF scheme of the homogeneous 
electron gas. This question will be investigated in the next extension of the work. 
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One point to be underlined is that in the homogeneous electron gas, it is known that the HF scheme predicts a 
non physical vanishing density of states at the Fermi level [44H48J , it is clear that this property is eliminated in there 
result of the iterative process in the proposed scheme. 

It should be also remarked that the present work, had been motivated, precisely by the idea of deriving in its further 
extension, the crystal symmetry breaking and spin-space entanglement effects obtained in Ref. 36] for La^CuO^. 
In these works, a simple model of the CuO planes was constructed. It was defined by a set Coulomb interacting 
electrons, for which the free Hamiltonian was taken as a Tight Binding one. The electrons were assumed to half fill 
a single Tight Binding band. The lattice of the model was chosen as the 2D crystal formed by the Cu atoms in the 
CuO planes. The Coulomb interaction in the system was then considered in the HF approximation. In the study, a 
phenomenological value of the dielectric constant of nearly e = 10 was assumed to screen the Coulomb interaction. 
This consideration was done in order to match the width of the HF spectrum of the Coulomb interacting electrons, 
with the width of the single half filled band crossing the Fermi level, arising in the Matheiss full band structure 
calculation in Ref. 29]. After that, by simply eliminating some usually imposed symmetry constraints on the HF 
procedure, a fully unrestricted HF solution of the HF problem predicted an insulator gap and the antiferromagnetic 
structure for the La2CuOi, which are considered as pure strong correlation effects in the material. 

The above circumstances, also motivated the idea of constructing a full HF band calculation scheme, eventually 
being able to predict the assumed dielectric screening in Ref. [36j . by also incorporating the mentioned crystal 
symmetry breaking and spin-spacial entanglement effects, in describing the Mott strong correlation properties of the 
LaiCuO^ material. The possibility to derive the above mentioned (just fitted in [36]) dielectric constant, is suggested 
by the following circumstance. The usual HF band structure evaluations of the La-iCuO^ are known to predict an 
enormous band gap of nearly 17 eV [37| . However, after dividing this energy gap by the value of 10 for the dielectric 
constant employed in the described model of CuO planes, gives a result of 1.7 eV, which is close to the measured gap 
of the LCL2CUO4 of 2 eV. Then, the expectation arises that a band calculation scheme based in the proposed modified 
HF scheme, could work in predicting the main strong correlation properties of the La^CuO^, and eventually of its 
close related transition metal oxides. The investigation of these possibilities is expected to be considered elsewhere. 
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